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' Abstract. We compute magnetic moments of baryons with a heavy quark in the bound state 

approach for heavy baryons. In this approach the heavy baryon is considered as a heavy meson 
bound to a light baryon. The latter is represented as a soliton excitation of light meson fields. We 
obtain the magnetic moments by sandwiching pertinent components of the electromagnetic current 
operator between the bound state wave-functions. We extract this current operator from the coupling 
to the photon field after extending the action to be gauge invariant. 
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^ : INTRODUCTION 

^ . This talk is dedicated to Joe Schechter on the occasion of his 65* birthday. Joe has 

^ [ made important contributions to the descriptions of hadrons with a heavy quark 

O I In particular he has successfully developed a hadron model that is valid for finite masses 

r-| ■ of the heavy quarks but naturally leads into the heavy quark effective model as these 

P-i! masses tend towards infinity LSfcBEl]. In this talk we will present an application of this 

p^l approach. 

The basic idea is to consider an effective meson model that includes both, light de- 
grees of freedom and mesons that contain a heavy quark. In this model the interaction 
of the light degrees of freedom is governed by chiral symmetry. On the other hand the 
^ ■ heavy quark components of the heavy meson field are coupled such that the Lagrangian 

. reproduces the heavy quark symmetry, i.e. spin and flavor independence of the interac- 

tions, as the masses of the heavy mesons tend to infinity. It is thus essential that this part 
of the model contains both the heavy pseudoscalar and heavy vector meson fields. 

It is well established that the light sector of such a model may contain soliton solutions 
to the classical equations of motion [6, /]. After appropriate quantization these solutions 
may be identified with baryon states as in the Skyrme model framework [81. Essentially 
this is a practical realization of the large Nc picture for baryons in QCD [9]. The heavy 
mesons interact with the light mesons according to the roles of chiral symmetry. Substi- 
tuting the soliton for the light meson fields into that interaction Lagrangian generates a 
background potential for the heavy meson field. This potential has bound state solutions 
that, when combined with the soliton, describe baryons with a heavy quark I IQ. 11..12ll . 
It is then, of course, interesting to study properties of such a heavy baryon by com- 
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puting matrix elements of pertinent operators. The most prominent operator that comes 
to mind is the one of the electromagnetic current. In particular the spatial components 
of this operator determine the magnetic moments of the considered baryon. We obtain 
this operator by incorporating the photon field such that the total action becomes gauge 
invariant. The current operator is then the coefficient that is linear in the photon field. 

This presentation is based on on-going studies regarding electromagnetic properties 
of baryons with a heavy quark in the above described bound state picture. These studies 
will be presented elsewhere lll3ll . Therefore the numerical results presented here should 
be considered preliminary. We will commence this talk with a brief a review of the 
soliton picture for the light mesons when vector mesons like p and O) are included. 
Such an extension has proven necessary in many aspects, as the proton-neutron 
mass difference lll4l] . t he p redicted pion-nucleon phase shifts lllSh or the axial singlet 
charge of the nucleon [16]. We will then discuss the occurance of the bound state and 
describe the quantization of the combined soliton bound state system that leads to the 
heavy baryon state. Finally we will compute the magnetic moments of heavy baryons as 
matrix elements of such bound states. 



THE MODEL LAGRANGIAN 

In this section we review the classical, i.e. leading order part in the 1 /Nc expansion of 
the bound state description for the heavy baryons in the soliton picture. 



Light Mesons 



For the sector of the model describing the light pseudoscalar and vector mesons we 
adopt the chirally invariant Lagrangian discussed in detail in the literature lfl7L l6ll. This 
Lagrangian can be decomposed into a regular parity part 
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tr [ivv (P)F^'' (P)] +mltT [R^R^] (1) 



and a part which contains the Levi-Civita tensor, e^vajS • The action for the latter is most 
conveniently displayed with the help of differential forms p = p^dx^, etc.: 
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In eqs dl]) and Q we have introduced the abbreviations 

p^ = i(^<9^^t_^t^^^)^ v^ = i(^^^^t + ^t^^^) and R^ = p^-\v^. (3) 



Here ^ refers to a square root of the chiral field, i.e. U = Furthermore F^v (p) = 
dfipv — dvpji — ig [pfi.pv] denotes the field tensor associated with the vector mesons p 

and CO, which are combined in = ^(O^I + p^T°j /2 when the reduction to two light 

flavors is made. The parameters ^,71, etc. can be determined (or at least constrained) 
from the study of decays of the light vector mesons such as p ^ 2;r or ft) ^ 3;r [6]. The 
complete determination of all parameters, however, also requires some information for 
the light baryon sector, see below. 

The action for the light degrees of freedom (/ ^ + Fg) contains static soliton solu- 
tions. The appropriate ansatze are 

^(r)=expf^r-TF(r)y ft)o(r) = ^ and p,,(r) = ^e,,«r, (4) 
V2 / g gr 

while all other field components vanish. The resulting non-linear Euler-Lagrange equa- 
tions for the radial functions F(r), ft)(r) and G(r) are solved numerically subject to the 
boundary conditions F{0) = —7t, ft)'(O) = and G(0) = —2 while all fields vanish at 
radial infinity |6]. These boundary conditions are needed to obtain a consistent baryon 
number one configuration. 



The Relativistic Model for the Heavy Mesons 



In this subsection we present the relativistic Lagrangian, which describes the coupling 
between the light and heavy mesons ||3|] 
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Here we have allowed the mass M of the heavy pseudoscalar meson P to differ from 
the mass M* of the heavy vector meson Q^. Note that the heavy meson fields are 
conventionally defined as row vectors in isospin space with P = (P^'Y ^D^P = (D^^P^y 
etc.. The co variant derivative introduces the additional parameter a: 



D^P^ 



{d^ - iagp^ -i{\-a)v^)P'^ = {d^- iv^ - igaR^) P^ , (6) 
D^Ql = {d^-iv^-igaR^)QX. (7) 
The covariant field tensor of the heavy vector meson is then defined as 



(8) 



The coupling constants d,c and a, which appear in the Lagrangian (|51), have still not 
been very accurately determined. In particular there is no direct experimental evidence 
for the value of a, which would be unity if a possible definition of light vector meson 
dominance for the electromagnetic form factors of the heavy mesons was adopted 
The other parameters in (O will be taken |2] to be: 

d = 0.53, c = 1.60; 

M = 1865MeV, M* = 2007MeV, D- meson; 

M = 5279MeV, M* = 5325MeV, B- meson. (9) 

It should be stressed that the assumption of infinitely large masses for the heavy 
mesons has not been made in ©. However, a model Lagrangian which was only required 
to exhibit the Lorentz and chiral invariances would be more general than the relativistic 
Lagrangian ©. Rather the coefficients of the various Lorentz and chirally invariant 
pieces in the relativistic Lagrangian dSJ) have precisely been arranged to become spin- 
flavor symmetric in the heavy quark limit M,M* — > oo 



Bound States 

Here we briefly review the origin of bound states in the S- and P-wave heavy meson 
channels. These orbital angular momentum quantum numbers refer to those of the 
pseudoscalar component (P^) of the heavy meson multiplet (P^, Qji). More details are 
provided in refs. B4j]. 

For the P-wave channel the appropriate ansatz reads 
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Note that here p refers to a properly normalized spinor which describes the isospin of the 
heavy meson multiplet. There exist also S-wave bound states {P^ = ^^4>(r)pe'^', . . .) 



/Atz 

which are, however, not relevant for the discussion of the magnetic moments. For more 
details see ref. fl]. 

We substitute the soliton background ^ for the light fields into the Lagrangian Q 
and find the equations of motion for the radial functions in the ansatz (flOb . Due to the 
soliton background there are solutions to these equations with |£| < M. These are the 
bound states we are looking for. 

As the relativistic Lagrangian ^ is bilinear in the heavy meson fields the resulting 
equations of motion are linear. Hence the overall magnitude of the solution is not fixed 
by the equation of motion. Nevertheless, the equations of motion for the heavy meson 
fields allow us to extract a metric for a scalar product between different bound states. 
In particular its diagonal elements serve to properly normalize the bound state wave- 
functions. The Lagrange function which results from substituting the ansdtze (El) and 



(ITUb may generally be written as 



L = -Mci [F, G, CO] +Ie [F, G, co; 4>, ^o, , ^2] P^(e)p(e) • (1 D 

Here Md denotes the soliton mass whose minimum determines the light meson pro- 
files F, G and CO. The explicit expressions for the functionals /g are given in ref. 0]. The 
subscript refers to the explicit dependence on the energy eigenvalues. Upon canonical 
quantization the Fourier amplitudes p(e) and p^(e) are respectively elevated to annihila- 
tion and creation operators for a heavy meson bound state with the energy eigenvalue £. 
Demanding that each occupation of the bound state adds the amount |e| to the total 
energy yields the normalization condition 
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in addition to the canonical commutation relation [p,(£), (e')] = dijd^ ^i. Note that for 
bound states the energy eigenvalues are discretized. For the P-wave channel we obtain 
the normalization condition 
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from eq (fT^ . For convenience we have employed the abbreviation i?^ = cosF — 1 -f- 
oc ( G — cosF) . The analogous condition for the bound state wave function in the S- 
wave channel is given in ref. [4]. The radial profiles associated with these normalizations 
are displayed in figure [l] for the choice a = —0.3. The parameters in the relativistic 
Lagrangian dSj) have been set to the charm sector Q. The parameters entering the light 
meson Lagrangian (1 1121) are given in eq (l25t . 



CRANKING THE BOUND HEAVY MESON STATE 

It can easily be verified that the field configurations for both the light mesons Q and the 
heavy mesons configurations (fTOI) obtained above are neither eigenfunctions of the spin- 
nor the isospin generators. Here we will construct such eigenstates from the soliton - 
bound state system. 



Collective Coordinates and Their Quantization 

In order to generate states which correspond to physical baryons a cranking procedure 
is employed. In the first step collective coordinates, which parameterize the (iso-) spin 
orientation of the meson configuration, are introduced via 

^-^A{t)^A'^{t) and ^ A{t)p^A'' (t) . (14) 



P - wave S - wave 

E = 1388 MeV E = 1618 MeV 
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FIGURE 1. The profile functions for the bound state wave-functions in the P-wave (left panel) and 
S-wave (right panel) channels. These functions are measured in units of my — 773MeV. See text for the 
specification of the remaining parameters. Figures adopted from ref. |4]. 



The time-dependence of the collective rotations is measured by angular velocities Q 

A\t)jA{t) = ^T-a. (15) 

In addition to the collective rotation of the soliton configuration, field components are 
induced that vanish classically. For the light vector mesons these are L2J 

COi = ^(p{r)£ijkajrk and = ^i{r)ak + ^2{r)r ■ ■ (16) 

The light meson Lagrangian then contains a term which is quadratic in the angular veloc- 
ities. The constant of proportionality defines the moment of inertia oc^ [F, G, (O; i^i , (^2, 9] • 
The radial functions (p{r), ^i{r) and ^lir) are obtained from a variational approach to 
10]. Here it is worth mentioning that is of order Nc- 

In analogy to eq (fT4ll the heavy meson fields also need to be rotated in isospin space, 
pt_^^(^)pt and Ql^A{t)Ql. (17) 

Substituting the collectively rotating configurations into the total Lagrangian finally 
yields (the subscript £ = 0, 1 refers to S and P wave bound states) 

Le = -M,i + I^J^p'^p + ^a^Q\^XeP^^-'^P ■ (18) 

(£) 

The quantity /g ' may be extracted from ref. [4]. It has already been employed to obtain 
the bound state profiles functions, e.g. eq. dTUt . The new quantity is the hyperfine 
parameter Xe whose explicit expression is also displayed in ref. fl]. 

Once the Lagrangian (fTSb for the coupling of the collective coordinates, A to the 
creation and annihilation operators, p of the bound is found, its quantization proceeds 



along the bound state approach to the Skyrmion In doing so, we have to 

construct Noether charges for spin and flavor. For this purpose we first have to consider 
the variation of the fields under infinitesimal symmetry transformations. For the isospin 
transformation we observe 
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Here refers to any of the iso-rotating meson fields and the ellipses represent terms, 
which are subleading in l/Nc, as e.g. time derivatives of the angular velocities which 
might arise from eq. (fT6t . Furthermore Dij{A) = (1/2) tr(T,ATyA^) denotes the adjoint 
representation of the collective rotations A. From eq (fT^ we conclude that the total 
isospin is related to the derivative of the Lagrange function with respect to the angular 
velocities 



Ii = -D,j{A) 



(20) 



Next we note that the total spin operator J contains the grand spin operator G 

Gi^Ji + Dr\A)Ij=Ji-jf . 
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with J" ' = dLi/dO.. As a consequence of the relation (l20b we have (7^°')^ = 7 = 
/(/+ 1). By construction the light meson fields do not contribute to G. Even more 
importantly and as has been noted before, the pieces of the heavy meson wave-functions 
(fTTl . which are located between the collective coordinates A and the spinor p, have zero 
grand spin too. With the normalization condition (fT3t one then finds 



(22) 



This relates the operator multiplying the hyperfine parameter in the collective La- 
grangian (fTSt to the spin and isospin operators. The collective piece of the Hamiltonian 
is obtained from the Legendre transformation 
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Here Lf^^ refers to the Q. dependent terms in eq dTSt . Finally the mass formula for a 
baryon with a single heavy quark becomes 



= Mel + + ^ 1) + (1 1)] 



(24) 



where contributions of ^{xj), which apparently are quartic in the heavy meson wave- 
function, have been omitted for consistency because terms of that order have been 
excluded from the very beginning. 



TABLE 1. The parameters of the mass formula i24\ for 
a =0.3. Here co = M — \e\ denotes the binding energy. 
The subscripts S and P refer to the bound states in the S 
and P wave channels, respectively. 

(Op Xp (Os Xs 

478MeV 0.114 247MeV 0.205 
713MeV 0.0045 573MeV 0.0052 



D-meson: 
B-meson; 



TABLE 2. Heavy baryons mass differences with respect to A^- or A^ for 
a — 0.3. Primes indicate negative parity baryons, i.e. S-wave bound states. All 
energies are in MeV. 
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From eq (EU we recognize that the spin degeneracy between baryons containing a 
heavy quark vanishes in the heavy quark limit because Xp approaches zero. Of course, 
this result is a direct consequence of the spin-flavor symmetry and would not have come 
out in case the various Lorentz and chirally invariant terms in eq ^ had been chosen 
arbitrarily. 

The parameters in the light sector, eqs. (1 1121) cannot completely be determined from 
properties of the corresponding mesons. The remaining (limited) parameter space is, 
however, more than fully constrained by a best fit to the mass differences of the low- 

lying ^ and | baryons in the light sector. This yields: 

^ = 5.57, mv=773MeV 
71 =0.3, 72 = 1.8, 73 = 1-2. (25) 

The resulting mass differences for the light baryons all agree within about 10% with 
experiment. The corresponding moment of inertia is a? = S.OOGeV ^ The parameters 
in the heavy baryon mass formula (l24b are listed in table [T] The resulting masses are 
displayed in table El We observe that the model represents a sensible picture of the 
spectrum of baryons with a heavy quark. This even more motivates the investigation of 
their properties, as e.g. the magnetic moments. 



ELECTRO-MAGNETIC CURRENT AND MAGNETIC MOMENTS 



In this part of the talk we would first like to describe the extraction of an electromagnetic 
current operator for bound state configuration. We then compute the magnetic moments 



by sandwiching that operator between the heavy baryon states, 



B{I,Iy,J,J,)=^ £ d/l 1 D'=':%{A)\^,G3) (26) 

that are obtained by coupling the collective coordinate wave-function (represented by 
the Wigner D-f unction, D{A)) together with bound state wave-function to total spin J 
according to eq. (ETT) . The prefactor ^ is a suitably chosen normalization constant. 

First we introduce the photon field £/fi in close analogy to refs. [6, 7] for the light 
sector and ref. ^ for the sector of the heavy mesons. The aim is to construct a gauge 
invariant Lagrangian =Sf ((^ , p , P, 2; £/) and to extract the current operator as the term that 
is linear in the photon field: 7^™- = d££' jds^^ I^m-o' ^^S^*- ^^'^tor, this exercise 

has been done in refs. [6, 7], while ref. |2] provides the electromagnetic generalization 
of the covariant derivatives (I6I7I) . Here we also need to gauge the derivatives in which 

can be done by writing I^ItI] P/x — ^ P/x +^^('^2'^^ — ^^Q^)^\i- Obviously this gauge 
term vanishes when light pseudoscalar fields are omitted ((^ = I). Therefore such a 
contribution did not arise in the work of ref. [2] as those authors were only interested in 
electromagnetic processes of heavy mesons that did not contain pions. To this end the 
electromagnetic current associated with the Lagrangian Q reads 
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+ie^ef,,p, {pCF-^{p)Q''^-Q'F^^{p)CP^) {11) 

where C = C + ^^{^Q^^ + t,^Qt,) and C is the charge of the heavy quark that is 
contained in the heavy meson (C = 2/3 and C = — 1/3 in the charm and bottom sectors, 
respectively). Substituting both the (rotating) soliton and the bound state wave functions 
we compute 

\ j ^e.m.) ^ = pisW^3 + Mko^33(^) + MsjP^ yP + Mki^ss (^)P V • (28) 

Here flsfl^ . . . , are functionals of all radial profiles that have been computed earlier. 
The functionals Uy n and /ivo do not contain the heavy meson wave functions and are 
listed in refs. |'7','l9]. Note that in the first term on the right hand side we made explicit 
the appearance of the moment of inertia, a? to make contact with earlier work [ 1 1. 12ll. 
The detailed form of the functionals }ls i and iiy i that contain the bound state wave 



functions is presented in the appendix, see also ref. 11131] . Using the quantization rules 
for the collective coordinates discussed above we can sandwich this operator between 
baryon states. 
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^B = :^{B\Jd'r(rx 7e.m. j ^ \B) (29) 



TABLE 3. Preliminary results for the magnetic moment of baryons with a heavy quark 
in nucleon magnetons. For comparison, we predict /i(n) = 2.58 and jj.{n) = —2.26 which 
compares reasonably well with the empirical data 2.79 and —1.91. 

A, E++ EO E- E+^A, a, EO EO E, E+ ^ A, 
0.12 2.56 0.17 -2.22 -2.39 -0.02 2.63 0.21 -2.19 2.41 



yielding lfTll[l2 1 

li{p) = 5^5,0 + fMv,o = 5Ms,o-fMv,o 

= 2A5,i 

) = 5^5,0 + fMv,o-^A5,i + fMv,i 

) = f M5,o - f Mv,o - ^A5,i - f Mv,i 
Ai(4=^^A2) = -§AiKO-fAiv,i (30) 

where jls^i = Xt^Sfi — t^s,\ arises from the quantization rule a^Q. = J^°^ — xG. In table|3l 
we present numerical results for the magnetic moments in eq. (OUb . Those results origi- 
nate from iisA = — 0.20(0. 05)n.m. and /iyj = — 0.04(— 0.01)n.m. in the charm (bottom) 
sectors. 

We note that, in particular for Ac, our results compare reasonably well to those 
predicted from a spectral sum rule analysis ['20]. On the other hand, our predictions 
for the transition magnetic moments Ln — > Ag are slightly larger (in magnitude) than 
those obtained from QCD sum rules MM- 

We furthermore observe that with increasing masses of the heavy mesons the param- 
eters M5.1 and i tend to zero. Then all magnetic moments are dominated by their 
light meson components f22]. Of course, this is merely a reflection of the heavy quark 
symmetry which implies that the magnetic moment of the heavy quark vanishes as its 
mass becomes infinite. 



CONCLUSIONS 

In this talk we have presented a study on the magnetic moments of baryons with a heavy 
quark. The picture for such heavy baryons is that of a heavy meson being bound to 
a baryon that is built out of light quarks. We adopt the soli ton picture for baryons in 
which light baryons are represented by soliton excitations of light meson fields. This is 
in the spirit of the Skyrme model but we have also included light vector meson degrees 
of freedom in the construction of the soliton. Such extensions by short range fields are 
needed for various phenomenological reasons. 

We have obtained the heavy meson bound state from a relativistic Lagrangian that 
does not manifestly reflect the heavy quark symmetry. Rather this model is developed 
for physical values of the heavy meson masses. However, in the (academic) limit of 
large masses it is demanded to reflect the heavy quark symmetry. This condition strongly 



constrains the number of free parameters. Also, it requires the model to contain heavy 
vector meson fields rather than just the pseudoscalar degrees of freedom as e.g. in the 
bound state approach to the Skyrme model and its application to the charm sector Jl^]. 
We have already observed in ref. that such extensions are also required to properly 
reproduce the spectrum of baryons with a heavy quark. 

Once the model has been set up, it is straightforward to extract the operator of the 
electromagnetic current by gauging the action with respect to the photon field. In turn 
we have then used the methods to of refs. ifTllfl?!] to compute the magnetic moments of 
baryons with a heavy quark. At this time the results are still somewhat preliminary and 
in future we wish to further explore the connection with the heavy quark limit for the 
magnetic moments [22J . 

We hope that Joe considers this study not only an interesting application of ideas that 
he decisively helped to develop, but also worth to be reported on the occasion of his 
birthday. 
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APPENDIX 

In this appendix we list the explicit expressions for the parameters in eq. 
the P-wave bound state wave-function (fTOb . 
For the isoscalar case we have 



that involve 
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ms{r) = -2{l + \Ra)^' 
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-t/fJf'»F2-sinF»FiWo + ^^ (rG)'»F2-2G%)<I>, (31) 
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where again Ra = cosF — 1 + a(l +G — cosF) and primes denote derivatives with 
respect to r. For the isovector contribution we find 
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+{ {%r + ^'2 + Ra'i'i) »Pi + (1 + ^/?«) - [/?«*Po + (e - f «)»P2] ^ojcosF 
+dsm2F^>o'¥i ~ {r(o'^2 - 2G'^>o) OcosF , 
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m\;'{r) = -2MdrsinF'¥2<i> 

+d [r (»Fo»P^ - %^2) + + 27?a»Fo»Fi + 2r(£ - ^)^i^2] (32) 

The analytic expressions for flsfi and flyfi can be found in refs. 
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